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Evolution of gaseous disk viscosity driven by supernova explosion 
in star-forming galaxies at high redshift 
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ABSTRACT 

Motivated by Genzel et al.'s observations of high-redshift star-forming galax- 
ies, containing clumpy and turbulent rings or disks, we build a set of equations 
describing the dynamical evolution of gaseous disks with inclusion of star for- 
mation and its feedback. Transport of angular momentum is due to "turbulent" 
viscosity induced by supernova explosions in the star formation region. Analyti- 
cal solutions of the equations are found for the initial cases of a gaseous ring and 
the integrated form for a gaseous disk, respectively. For a ring with enough low 
viscosity, it evolves in a slow processes of gaseous diffusion and star formation 
near the initial radius. For a high viscosity, the ring rapidly diffuses in the early 
phase. The diffusion drives the ring into a region with a low viscosity and start 
the second phase undergoing pile-up of gas at a radius following the decreased 
viscosity torque. The third is a sharply deceasing phase because of star formation 
consumption of gas and efficient transportation of gas inward forming a stellar 
disk. We apply the model to two z ~ 2 galaxies BX 482 and BzK 6004, and find 
that they are undergoing a decline in their star formation activity. 

Subject headings: galaxies: evolution - galaxies: high-redshift 



1. Introduction 

Exciting progress has been made in observing assembly of gas and star formation in 
galaxies from low to high redshifts. Optically/UV-selected star-forming galaxies at z ~ 
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0.5 2.5 exhibit a fairly tight relationship between stellar mass and star formation rate, 

which is strongly suggestive of star formation happening with a high duty cycle (Noeske et 
al. 2007; Elbaz et al. 2007; Daddi et al. 2007). This may relate with random accretion 
onto supermassive black holes (SMBHs) triggered through minor mergers located at galactic 
centers (Wang et al 2009), showing episodic activities of the SMBHs (Wang et al. 2006, 

2008) . Major mergers are undoubtedly taking place in very luminous submillimeter galaxies 

at z ~ 1 3 (Tacconi et al. 2006; 2008), but only a small fraction of this sample (~ 1/3) 

is undergoing obvious major mergers from the images of the sample (Foster Schreiber et al. 
2006; Genzel et al. 2006; Wright 2007; Law et al. 2007). Detail studies of Spectroscopic 
Imaging survey in the Near-infrared with SINFONI (SINS) of high-redshift galaxies show 
that about 1/3 of the samples are rotation-dominated turbulent star-forming rings/disks; 
another one-third are compact and velocity dispersion-dominated objects and the remains 
are clear interacting/merging systems (Genzel et al. 2008; Foster Schreiber et al. 2009). 
These dynamical evidences show that rapid and continuous gas accretion via "cold flows" 
and/or minor mergers are playing an important role in driving star formation and assembly 
of star- forming galaxies at 2; > 1 (Genzel et al. 2008). Interestingly, recent numerical 
simulations produce results agreeing with the observations (Bower et al. 2006; Kitzbichler & 
White 2007; Naab et al. 2007; Guo & White 2008; Dave 2008; Brooks et al. 2009; Dekel et al. 

2009) . Minor mergers as a supply of gas may trigger star formation and subsequently ignite 
the SMBHs. The role of supernova explosions (SNexp) in driving interstellar turbulence 
and the existence of a supernova rate-velocity dispersion relation has been highlighted using 
numerical simulations (e.g., Wada & Norman 2002; Dib et al. 2006), and observationally by 
the spatial correlations of the inclinations of molecular clouds with respect to the Galactic 
disk attributed to the effects of SNexp (Dib et al. 2009; see an extensive review of Mac Low 
& Klessen 2004), and by the existence of a starburst-AGN connection (Chen et al. 2009). 
SNexp is potentially playing a role in the dynamics in star-forming galaxies. 

Distributions of density and angular momentum of a gaseous disk are generally con- 
trolled by the viscosity, if driven by viscosity of turbulence excited by SNexp, stellar distri- 
bution will be determined by this feedback process (Lin & Pringle 1987). Steady analytical 
solution of circumnuclear disk of the Galaxy with injection of SNexp energy has been ob- 
tained by VoUmer & Beckert (2003), coevolution of SMBHs and their hosts could be driven 
by SNexp-excited turbulence (Kawakatu & Wada 2008) as well as in the starburst-AGN con- 
nection in light of strong correlation between the Eddington ratios and speciflc star formation 
rates (Chen et al. 2009). Figure 1 shows a clear trend of concentration of galaxies with star 
formation rate densities in the SINS sample given by Genzel et al. (2008), strengthening the 
interplay and regulation between SN explosions, cold gas accretion, and minor mergers. 

In this Letter, we consider the secular evolution of gaseous disk, which could be formed 
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during minor mergers, by including feedback from SNexp in star formation or burst regions 
characterized by the Ha ring in the SINS sample. We find the analytical solutions, which 
show two kinds of properties: 1) piled-up gas and 2) fast diffusion depending on the initial 
conditions of gas density. We apply the models to high-redshift galaxies. 
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Fig. 1. — Mass concentration of galaxies defined by a ratio of dynamical mass within 0".4-l".2, Cm = 
-^dyn(< 0".4)/Mdyn(< l"-2), Correlates with star formation rates for the SINS sample given in Table 2 in 
Genzel et al. (2008). This correlation strengthens the role of SNexp. 



2. Dynamical equations 



Considering the mass dropout of gaseous disk due to ongoing star formation, the mass 
conservation equation reads 

d 



(1) 



where Vr is the radial velocity of the gas, S is the surface density of gas, and is the 
surface density of star formation rates. We neglect mass injection of stellar winds back to 
the gaseous disk or escapes of winds from galaxies. Star formation process removes some 
angular momentum of gas, thus the conservation equation of angular momentum is given by 



(2) 



where Q = 2ttR^i'T, (dfl/dR) is the viscosity torque, u is kinematic viscosity, and fl is the 
angular velocity. Combining eqs. (1) and (2), we have 
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In this Letter, we explore the role of viscosity driven by SNexp in the secular evolution of 
gaseous disk. Similar to Kawakatu & Wada (2008), we assume that the vertical structure is 



-4- 



supported by turbulence pressure 

Ptur ( GM.H , \ 

^ = Pg(^^ + -GEj^p,-, (4) 

where M, is the black hole mass, G is the gravity constant, Q, — \f2a/R and ttGE = cr^/i? 
are used (eq. 1 in Thompson et al. 2005), Ptur = Pg^tur is the turbulence pressure, a is 
the dispersion velocity of galaxies determined by the total potential of gas, stars, and dark 
matter halo, and H is the half-thickness of the disk. The validity of solutions presented here 
is limited by the potential used in eq. (4), namely, GM,H/R^ < c'^/R, we thus have, R > 
GM,{H/R)/a^ ~ 10 Ms(r^oo{H/R)pc, where Mg = M./IO^Mq and CT200 = (r/200 km s-\ 
The energy equation of turbulence excited by SNexp is given by Wada & Norman (2002) 

Pg^ur Pg^ur /-C IT" 

—7 — = — Tr~ ^ ^'^★-^SN, [o) 

where tdis = -f^/Kur is the timescale of the turbulence, ^ is the efficiency converting kinetic 
energy of SNexp into turbulence, S*^ is the SNexp rate, and ii^sN is the SNexp energy. SNexp 
rate strongly depends on the initial mass functions, but we absorb its uncertainties into the 
parameter ^ and simply let be the star formation rate. 

The viscosity law z/ = aVturH is assumed commonly, where a is the viscosity parameter 
(Shakura & Sunyaev 1973). The Kennicutt-Schmidt's law displays oc S^'^, however, there 
are growing evidence for a linear relation for denser environment in star-forming galaxies from 
numerical simulations (Dobbs & Pringle 2009; Krumholz et al. 2009). We use = c^S 
for star- forming galaxies at high-redshift, where is a constant. We thus have the energy 
equation T/^^^. = ^c*if E'sn, where is the efficiency converting gas into stars. With the help 
of equations (4) and (5), we have 

u{R) = CP^ = 0.128 Co.usRt kpc' Gyr-\ (6) 

where ( = a {^c^Es^^f /cr^ = 0.128ao.i (C-sC-s-^si)^ (T^oo kpc'^Gyr'^ Ri = R/1 kpc, ao.i = 
a/0.1, = EsN/lO^^erg, = ^/IO'^Mq^, and = cj3 x 10"^ yr^^ Here, the star 
formation efficiency is assumed to be a constant independent of the gas density and radius. 
Such strong dependence of viscosity on the radius means that transportation of angular 
momentum will be much more efficient at outer radii than at the inner. 

The star surface density can be obtained from 

E,(P,i)= [ t^dt'^ [ cj:{R,t')dt' . (7) 



Vertical self-gravity may directly control the star formation as suggested by Thompson 
et al. (2005). We test self-gravity of the disk through Toomre parameter Q, defined by 
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Q = CsKq/tcGH oc where = 4fi^ + dQ'^/dlnr is the epicychc frequency, Cg = 

1.0X2^''^ km s~^is the sound speed, and T2 = T/IO^K is the gas temperature (Kawakatu & 
Wada 2008). We should keep in mind that there are two different aspects from the classical 
gaseous disk investigated by Lynden-Bell & Pringle (1974): (1) radius-dependent viscosity 
driven by SNexp and (2) dropout of mass and angular momentum due to star formation on 
the disk. These drive some interesting properties of gaseous disks. 

3. Solutions: evolutionary properties 
3.1. Gaseous ring 

A realistic initial condition is a single gaseous ring, which can be formed via one minor 
merger (Hernquist & Mihos 1995). Equations (3-6) have analytical solutions for the case of 
a single ring. After some algebraic operations (e.g., Kato et al. 1998), we have the analytical 
solution for the case 
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Fig. 2. — Properties of the surface density I](r, r) as to qq. Here, rp and Tp are solutions of dT./dr = and 
d'E/dr = 0. S(rp,Tp) is the peak or valley values at rp and Tp. The red and blue points correspond to peak 
and valley of the density, respectively. The stars represent rp. When qq > = 0.544, there is no solution, 
namely, the density monotonously decreases with time and radius. All the hollow symbols are the valley of 
the envelope of the E-evolution track whereas the solid ones are its peak. 



where Sq = Mq/ttRI, t = c^^t, r = R/Rq, and Mq is the total mass of the initial ring at 
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a radius Ro, and the dimensionless parameter qo is defined by 

?o = ^ = 0.584 a^lcZlRo'i (^-3^51)"' a^,,, (9) 

where i?o,i = -Ro/10 kpc. 

Comphcated evolutionary behaviors of the surface density are determined by go values as 
shown in Figure 2. For a given < qc = 0.544, there are two sets of roots of 9E(r, T)/dr = 
and (9S(r, T)/dr = 0, corresponding to the peak and valley of the envelope 

of S-evolution track at rp and Tp. Otherwise there is no root of the two equations, 
namely there is no peak or valley of the envelope. Evolution of S(r, r) and S*(r, r) can be 
found in Figure 3(a) and (b), respectively. We show three types of the solutions: (1) go < Qc 
(2) go = gc, and (3) go > gc- The go > gc solutions are quite simple and similar to the pure 
diffusion of a gaseous ring described in Lynden-Bell & Pringle (1974). The ring dramatically 
spreads over the space and forms stars. A stellar ring forms then as shown in Figure 3(b) 
go = 2 panel. For go = gc solution, there is a flat envelope without peak and valley, forming 
a very broader stellar ring as shown in the Figure 3(b) go = 0.544 panel. 

The go < gc solutions have complicated behaviors as shown in Figure 3(a) go = 0.1 
panel. Evolution of the gas ring can be divided into three phases. At early time, star 
formation is only important in the role of transportation of angular momentum, and the 
mass dropout can be neglected. Behaviors of the ring are very similar to the pure fluid ring. 
The gas rapidly diffuses inward mostly and outward a little bit carrying away the angular 
momentum. Though the star formation rates are very high, the total formed stars are not 
many since the duration of this phase is quite short. The second phase begins when mass 
dropout due to star formation becomes important. Decrease of the gas density is leading 
to decrease of star formation and weakening the viscosity torque of SNexp, resulting in 
accumulation of gas at the radius rp. The third phase starts when the viscosity torque is 
enhanced in light of the gas accumulation. This causes intensive star formation and exhausts 
gas, giving rise to an exponential decreases of gas density. 
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Fig. 3. — Evolution of gaseous ring and star surface densities for different parameter gg. In top and 
middle panel, black lines show the evolution of gas and star disk from ti = 0.005 to T2 = 5 with an interval 
AlogT ~ 0.15. Red, green, blue, cyan, and pink lines in bottom panel represent the Toomre parameter Q 
at r = 0.005,0.03,0.2, l,anc?5, respectively, which correspond to the same meanings in panels (a) and (b). 
We set the surface density of gas is So — IOOM0 pc~^ to calculate the Toomre parameter for three different 
go values. 

To understand such behaviors of solutions, we compare the timescales of star formation 
(t* = 1/c*) and gas advection (tadv = R/Vr ~ R^/v = 1/C,R?). According to the go- 
definition, it follows ta_dv/t* = 4go- When t^, <^ tadv at Rq (roughly corresponding to go 3> Qc), 
most of the gas is converted into stars at Rq and only a little gas is conveyed inward. In 
contrast, in the S> tadv case (namely, go -C gc) at initial radius Rq, star formation is so 
slow that most of the gas is advected inward. For a case with initial go <^ gc, most of gas 
is advected inward until tadv > 4- Pile- up of gas appears at Vp then and forms stars. The 
complicated behaviors of gas density are results as a consequence of competition between 
the two processes. 

Figure 3(b) shows the surface density of stars for the corresponding gas density. For 
go = 0.1, the gas rapidly diffuses inward and there is no significant pile-up of stars at the 
initial radius Rq. A stellar ring forms at roughly ~ 0.15i?o far away from its initial radius. 
The ring becomes broader with increases of go. When go = gc, there begins to be a tiny 
pile- up of stars at the initial radius. For a case with a large go, pile-up of stars happens very 
close to the initial radius and it becomes more conspicuous with go. The stellar rings become 
wider with go. Figure 3(c) displays evolution of the Toomre parameter Q, clearly showing 
Q < 1 in the star-forming region. Q-values are quite different in different regions. This 
strengthens the importance of self-gravity in the secular evolution of the gaseous disk. Since 
Q oc S~ Q basically follows the evolution of Egas- A generic property of the Q-parameter is 
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that it increases with time since the gas is being converted into stars. This is totally different 
from the steady gaseous disk with star formation under the presumed condition of Q = 1 
suggested by Thompson et al. (2005). 
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Fig. 4. — Evolution of gaseous disk and star surface densities for different initial disks. In the top and 
middle panel, black lines show the evolution of gas and star disk from ri = 0.005 to T2 = 5 with an interval 
Alogr = 0.15. Red, green, blue, cyan, and pink lines in bottom panel represent the Toomre parameter Q 
at T = 0.005, 0.03, 0.2, 1, and5, respectively, which correspond to the same meanings with panels (a) and (b). 
We set the surface density of gas is Eq = IOOMq pc~^ to calculate the Toomre parameter for three different 
7 values. The yellow lines in (a) are the initial surface densities. 

We would like to point out that the feedback from SNexp will be enhanced, if we use 
the non-linear Kennicutt-Schmidt's law to calculate the star formation rate, leading to more 
efficient viscosity to transport angular momentum. However, the main features given by eq. 
(8) will remain. A detailed discussion on this nonlinear effects is beyond the scope of this 
Letter, but it would be interesting to solve the nonlinear equations to display the effects in 
future work. This model may apply to different kinds of galaxies. 

3.2. Gaseous disk 

Successive multiple minor mergers may provide an initial gas widely distributed in the 
major galaxies (Hernquist & Mihos 1995; Bournaud et al. 2007). For a simplicity, we 
assume an initial distribution with a power law of radius as 0) = Sq {R/RaY, where 
R^ = icjACf' = 7.65 c_/^Co.u! kpc, So is the surface density at R^, and 7 is the index. 
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We integrate the Green's function over the radius range of the disk, and have 

^^r.,r).r^^(^-.-'S).--ii^-iy,r,, (10) 

Jo (7rr)2 J 

where r„ = R/Ra and rg = Ro/Ra- 

Figure 4 shows the evolution of gaseous disk. We calculate the cases of 7 = —1.0, 0, and2.0 
between Ri = 0.2Ra and R2 = 2.0R„. For 7 = —1.0 disk, the most dense part of the disk 
is located in the radius with < tadv, so the gas is converted to stars locally. In contrast, 
for 7 = 2.0 case gas is concentrated outside the region with > tadv The gas diffuses quite 
rapidly and transfers to the inner radius and forms stars, leaving a relatively broader stellar 
ring or disk. An initial homogeneous disk has a complex behaviors between the two cases 
above, inner part of the gas forms stars locally whereas outer part flows inward and then 
forms stars. The Toomre parameter shows strong variability with time, but holding Q < 1 
in most star-forming region. 




fl{kpc) fl(kpc) 

Fig. 5. — All the red thick hnes are the current status of galaxy BX 482 and BzK 6004. Given a = 
210km s^^ from observations for the model parameters of BX 482, we relatively arbitrarily choose Rq — 
lOkpc, then make out qq — 0.86. We have Sq = 9.3Mq pc~^ from star formation rates. For Bzk 6004, we 
find that the initial parameters are So = GQMq pc^^, Ri = 1.6, i?2 = 14.2kpc, 7 = 2, and t — 3.3 x lO^yr. 
The solid lines represent the past status of the two galaxies and dotted lines do there future status. 

The instability of the gaseous disks can be simply justified by the operating mechanism 
of the viscosity. When gas is piled up, the higher star formation rates can efficiently remove 
the angular momentum of gas, and then suppress the star formation rates. The system 
should be self-organized and stable. However, we should give a detail description of stability 
analysis in a future paper. 
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4. Applications 

Dynamics of the gaseous disk presented here allows us to compare current model with 
Genzel et al.'s observations of star- forming galaxies at high-redshift. Observations of SIN- 
FONI on ESO VLT provide details of the "cold gas" accretion. We choose two galaxies as 
illustrations to apply the present theoretical model. BX 482 is a ring-dominated star-forming 
galaxy, of which the ring radius is Rrmg — 7.0 ± O.Skpc with a narrower width AR ~ 1 kpc, 
its dispersion velocity is cr = 210km and the surface density of star formation rates is 
0.57 ± O.SMq yr~^ kpc^^. We hence choose the ring model (eq. 8). The bright Ha ring 
limits the lifetime of the ring, which is about ~ lO^yr in light of OB stars. The left column 
of Figure 5 gives the dynamical evolution of the ring. The Toomre parameter Q shows that 
the ring is self-gravity dominated. A stellar ring formed in the gaseous ring is still growing 
shown by the red line, but the Ha ring is becoming dimmer. 

BzK 6004 is also a rotation-dominated galaxy with relatively faint but broader Ha ring 
{AR/R ~ 1.0 from Figure 5 in Genzel et al. 2008). The ring radius is Rriug = 6.9±0.8kpc and 
the dispersion velocity of galaxy is o" = 240km s~^ and the surface density of star formation 
rate is 0.6 ± O.SM© yr^^ kpc~^. The current data does not allow us to distinguish the case 
of different 7 we thus just use 7 = 2 as an illustration. The dynamical evolution of the 
galaxy is shown by the red line in the right column of Fig. 5. The current stage of the ring 
is approaching the end states of the star forming. The stellar ring is still growing fast and 
may remain as a wide stellar ring. The Q- value shows that self-gravity is also dominated. 

We note that the solutions presented here tend to show an exponential stellar disk 
if the time is long enough. This agrees with the steady disks known for quite long time 
from Lin & Pringle (1987). The simple apphcation of the current model shows here that 
SNexp does indeed play a key role in the dynamical evolution of the gaseous rings or disks. 
However, the application should be improved with the more detail data. For example, the 
good-quality data of Ha ring profile and the spatially resolved spectroscopy will allow us to 
get the evolution stage of gas and star formation from the stellar synthesis. Particularly, the 
very sharp inner edge of the ring predicted here is caused by SNexp. This reflects feedback 
of the star formation. Additionally, application of the present model to a large sample (e.g. 
Foster Schreiber et al. 2009) will produce the initial conditions of the sample for a statistics 
in future. This will be invaluable to justify the origin of the cold gas. 

We would like to stress here that the present applications focus on the main features of 
galaxies at high-redshift. The ring or disk models for the two galaxies are less conclusive. 
We need more information of the galaxies, such as, stellar ages at different ring radii, so as 
to give more robust pictures. Degeneracies of parameters limit applications of the model. 
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5. Conclusions 

We set up a model of gaseous disk with star formation, in which the viscosity is driven 
by turbulence excited by SNexp. We get the analytical solutions for the initial condition of 
single gas ring, showing complicated secular evolutionary behaviors. Several different kinds 
of stellar distributions can be formed, including exponential types. The present model is 
able to reproduce starburst rings in two galaxies as illustrations. 

The model will be improved in several aspects in future by using: (1) nonlinear star 
formation law, (2) separated equations of stars and gas, and (3) initially twisted gaseous 
disk (Lu & Cheng 1991). These will help us to get the initial conditions as a consequence of 
minor mergers and allow us to study the evolution of stellar disks and bulges. 

The referee is greatly thanked for a very helpful report improving the manuscript. Prof. 
J.-F. Lu is acknowledged for useful discussions. We appreciate the stimulating discussions 
among the members of IHEP AGN group. The research is supported by NSFC-10733010 
and 10821061, CAS-KJCX2-YW-T03, and 973 project (2009CB824800). 
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